Spin-dependent WIMP scattering off nuclei 
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Chiral effective field theory (EFT) provides a systematic expansion for the coupling of WIMPs to 
nucleons at the momentum transfers relevant to direct cold dark matter detection. We derive the 
currents for spin-dependent WIMP scattering off nuclei at the one-body level and include the leading 
long-range two-body currents, which are predicted in chiral EFT. As an application, we calculate the 
structure factor for spin-dependent WIMP scattering off 129,131 Xe nuclei, using nuclear interactions 
that have been developed to study nuclear structure and double-beta decays in this region. We 
provide theoretical error bands due to the nuclear uncertainties of WIMP currents in nuclei. 

PACS numbers: 95.35.+d, 12.39.Fe, 21.60.Cs 



Introduction.- Cosmological and astrophysical obser- 
vations have established that more than 20% of the en- 
ergy density of our universe is dark matter, a weakly- 
interacting nonbaryonic form of matter, whose specific 
composition remains unknown [l[ . Promising candidates 
are weakly interacting massive particles (WIMPs), such 
as neutralinos, the lightest supersymmetric particles pre- 
dicted by extensions of the standard model. This has 
spurred direct detection of cold dark matter via elas- 
tic scattering off nuclei, requiring detailed knowledge 
of the response to WIMP induced currents in nuclei. 
This presents a challenging problem, because even if the 
coupling of neutralinos (or other particles) to quarks is 
known, it needs to be evaluated at the nucleus level in the 
nonperturbative regime of quantum chromodynamics 2] . 

Chiral EFT provides a systematic expansion for nu- 
clear forces and the coupling to external probes for mo- 
menta of order of the pion mass, p ~ ~ 100 MeV, 
which are typical momentum transfers in direct dark 
matter detection. In this Letter, we focus on spin- 
dependent neutralino scattering off nuclei, which is used 
to constrain WIMP properties |3| and is particularly sen- 
sitive to nuclear structure. We derive the currents for 
spin-dependent WIMP scattering at the one-body level 
and include the leading long-range two-nucleon currents, 
which are predicted in chiral EFT. The same physics 
is key for providing accurate predictions of neutrino 
breakup of the deuteron for SNO 0] and for Gamow- 
Teller and double-beta decays 

We apply the developed chiral EFT currents to calcu- 
late the structure factor for spin-dependent WIMP scat- 
tering off 129 ' 131 Xe, as xenon isotopes provide the tight- 
est limits on WIMP couplings @. To describe these nu- 
clei, we employ the largest many-body spaces accessible 
with nuclear interactions previously used to study nuclear 
structure and double-beta decays in this region Q. 

Spin- dependent interactions. -We focus on spin-depen- 
dent WIMP scattering, for which the low-momentum- 
transfer Lagrangian density C is taken to be a spatial 



axial- vector-axial- vector coupling 0, Q , 
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where Gp is the Fermi coupling constant, x the neu- 
tralino field, ip q the fields of g = u,d,s quarks, and A q are 
neutralino-quark coupling constants. At the one-nucleon 
level, the quark currents are replaced by their expecta- 
tion value in a single nucleoli. This leads to one-body 
(lb) axial-vector currents Jij_b. Summing over all A nu- 
cleons in a nucleus, one has 
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where p = pj — p' t denotes the momentum transfer from 
nucleons (with mass m) to neutralinos. The isoscalar 
part J° lb with coupling a — (A u + Ad)(Au + Ad) + 
2 A s As receives contributions from the isoscalar combi- 
nation of the u and d quarks to the spin of the nu- 
cleoli, as well as from the s quark. Here Ag denotes 
the matrix element of 4> q 1l5'4>q i n the nucleon up to the 
spin a/2 The isovector coupling can be written as 
ai = (A u — Ad)(Au — Ad) — (A u — Ad)gA, with the axial 
coupling constant gA- Therefore the isovector part J 3 lb 
of the axial-vector WIMP coupling to the nucleon is up 
to a factor i?^ identical to the axial-vector part of the 
weak neutral current. 

Chiral EFT and WIMP currents.- The weak neutral 
current was derived within chiral EFT for calculations of 
low-energy electroweak reactions. At order Q° and Q 2 , 
there are only one-body currents, which for the isovector 
part of the axial- vector WIMP current lead to 
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In chiral EFT, the p dependence of the axial and pseudo- 
scalar couplings gA(p 2 ) and gp(p 2 ) is due to loop correc- 
tions and pion propagators, to order Q 2 [Toj : 
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with A A = 1040 Me V, pion mass = 138.04 MeV, 
pion decay constant F n = 92.4 MeV, and g vpn — 13.05. 
In previous calculations of WIMP scattering off nuclei, 
the smaller 1/A 2 ^ terms are generally neglected and the 
Goldberger-Treiman relation is implicitly used to write 

f£ « ^ (compare with Ref. 1). 

Neglecting the strange quark contribution, the neu- 
tral axial-vector currents are isovector in the Standard 
Model. Therefore, knowledge of the isoscalar currents 
in the nucleon is based on models. These suggest that 
the Q 2 one-body currents have a form-factor mass-scale 
~ A a [ll[ . Because the isovector 1 / A\ terms contribute 
at the few percent level to spin-dependent WIMP scat- 
tering, we do not attempt to model these and neglect 
higher-order isoscalar current contributions. 

Two-body currents.- At order Q 3 , two-body (2b) cur- 



rents enter in chiral EFT [12j. For spin-independent 
WIMP scattering, the potential impact of meson- 
exchange currents has been pointed out in Ref. [13|. In 
chiral EFT, the long-range parts of 2b currents are pre- 
dicted. Due to their pion-exchange nature, they only 
have an isovector part, J 2 b = J2t<j^ij- For the axial- 
vector weak neutral current, this leads to 
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where t 3 = (j\ x t 2 ) 3 and the same for er x , k = ^(p' 2 — 
P2 - Pi + Pi) and q = ±(pi + pi - p 2 - p 2 ). The low- 
energy couplings C3 , C4 relate different processes and are 
determined in the pion- nucleon or NN systems 14 1. 

Following Ref. |5|, we include the normal-ordered lb 
parts of chiral 2b currents. This is expected to be a very 
good approximation in medium-mass and heavy nuclei, 
because of phase-space arguments [H| . We sum the sec- 
ond nucleon over occupied states in a spin and isospin 
symmetric reference state or core: J° 2b = Ylj0-~ Pij)^%^ 
where is the exchange operator. Taking a Fermi gas 
approximation for the core we obtain the normal-ordered 
lb current. This contains two terms. First, it leads to 
an in- medium renormalization of the axial coupling 
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Fermi momentum, and I{p, P) is due to the summation 
in the exchange term, 
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Such a renormalization is also expected considering chi- 
ral 3N forces as density-dependent two-body interac- 
tions 16], We take a typical range for the densities 
in nuclei p = 0.10.. .0.12 fm" 3 . This leads to I(p,P) = 
0.64. ..0.66, using the Fermi-gas mean-value P 2 = 6fcf/5, 
and I(p,P) = 0.58. ..0.60 for P = 0. So the total- 
momentum dependence is very weak, and we therefore 
take P — in the following. The resulting 2b contribu- 
tion to the axial-vector WIMP current can be included 
as a density-dependent renormalization a±(l +Sai), with 



^ = -A 7(p , P = 0) A (2c4 _ C3) + _L 
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Normal ordering leads to a second contribution which 
renormalizes the pseudo-scalar coupling. At P = 0, this 
takes the form for the weak neutral current 
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For neutrinoless double-beta decay, we approximated (p • 
at) p — > p 2 /3(Ti, which neglects the small tensor part 
between + ground states [H . For the axial- vector WIMP 
current, we fully include the contribution of Eq. ([9]) via 
a density- and momentum-dependent modification of the 
pseudo-scalar part. To this end, we define 



6a? {p 2 ) = -2c 3 



P 



F 2 4m 2 + p 2 



(10) 



In addition to the long-range pion-exchange currents, 
there are short-range 2b currents both for the isoscalar 
and isovector parts, which are included as contact terms 
in chiral EFT. For Gamow- Teller transitions, which are 
an isospin rotation of the axial-vector WIMP currents, 
the contributions from long-range 2b currents were found 
to dominate over the short-range parts in medium-mass 
to heavy nuclei (for typical short-range couplings cd), 
because c 3 ,C4 are large in chiral EFT without explicit 
Deltas [f|. Therefore, we neglect the short-range 2b 
current contributions here, consistently with neglecting 
higher-order lb isoscalar currents, as discussed above. 

WIMP-nucleus scattering.- The differential cross- 
section for spin-dependent WIMP elastic scattering on 
a nucleus in the ground state with total angular momen- 
tum J is given by Q 



with total momentum P = p; + p^. Here, p = 2fcp/(37r 2 ) 
is the density of the reference state, /cf the corresponding 
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where v is the WIMP velocity and Sa(p) is the axial- 
vector structure factor, which can be decomposed as 

Mp) = E (\(J\\TlHp)\\J)\ 2 + \(J\\£1(p)\\J)\ 

Lodd 

. (12) 

The sum is over multipoles L with reduced matrix ele- 
ments of the transverse electric T£ l 5 and longitudinal C\ 
projections of the axial- vector currents. These can be ex- 
panded in vector spherical harmonics [|, 17 1, leading to 
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Expressions for M L ^(pv.j) = j l> (pri)[Y L > (i-i) <Ti\ h (with 
V and cr coupled to L) are given in Ref. [13] • The pseudo- 
scalar currents do not contribute to the transverse elec- 
tric multipoles, and for the longitudinal part, one can 
replace (p • <Ti) p — > p 2 <Ji. As a result, the p 2 /A A terms 
in the longitudinal response cancel [see Eq. (d])] and only 
contribute to the transverse electric multipoles. 

Nuclear structure and results.- We carry out state-of- 
the-art large-scale shell-model calculations for the struc- 
ture of 129 Xe and 131 Xe. The valence space for both 
protons and neutrons comprises the Ogi/2, la's/ 2, W3/2, 
2s!/ 2 , and 0/in/ 2 orbitals on top of a 100 Sn core. In the 
case of 129 Xe, the number of particles in the ld 3 / 2 , 2,s 1/ / 2 , 
and 0/in/2 orbitals was limited to 3, in order to make 
the calculations feasible (the matrix dimension for this 
space is 3.5 x 10 8 ). We have used the so-called gcn5082 
interaction Q, which is based on a G- matrix with em- 
pirical adjustments, mainly in the monopole part, to de- 
scribe nuclei within this region. The same interaction 
and valence space have been used to study structure and 
double-beta decays. Calculations have been performed 
with the shell-model code ANTOINE 



18] 



Figure Q] shows the resulting spectra for 129 Xe and 
131 Xe, measured from the ground-state energy. In both 
cases, the experimental ground state and the overall or- 
dering of the excited states are remarkably well described, 
which represents a clear improvement with respect to pre- 
vious work This reflects the quality of the inter- 
action and valence space used. The resulting magnetic 
moments, with standard g factors for this region 20j, are 
fi = — 0.72^/v and 0.86^/v for 129 Xe and 131 Xe, respec- 
tively, which agree within ~ 20% with experiment. How- 
ever, we emphasize that the magnetic moments probe 
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FIG. 1: (Color online) Comparison of calculated spectra of 
129 Xe and 131 Xe with experiment. The interaction and va- 
lence space are the same as in structure and double-beta de- 
cay studies in this region 0. The 131 Xe diagonalization is in 
the full space, while the 129 Xe case was restricted (see text). 



different physics, as the axial- vector WIMP couplings at 
p = are an isospin rotation of the Gamow- Teller op- 
erator, not a coupling to magnetic moments. In fact, 
one could fine-tune the effective g factors to improve the 
agreement with magnetic moments (which is sometimes 
pursued [19j), but this leaves the predicted axial- vector 
structure factor unchanged. 

In the limit of low momentum-transfer, p — 0, the 
axial-vector structure factor reduces to the expectation 
values of the total proton and neutron spin operators, 

^(°) = -^\(ao + a\)(J\\S p \\J) + (a - a\)(J\\S n \\J)\ 2 , 

(15) 

where a[ = a\(l + 5a\) includes the effects from chiral 
2b currents. For the expectation values (S n , p ) = (JM = 
J|S 3 p | JM = J), we obtain 
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(16) 



As expected for odd-mass nuclei with even number of pro- 
tons, (S„) 3> (S p ). The predicted spin expectation values 
are qualitatively similar to those of Ressell and Dean [2l[ , 
who considered two different (Bonn A/Nijmegen II) in- 
teractions, although their (S p ) = 0.013/0.028 for 129 Xe 
is larger and their (S„) values for 131 Xe are 20% smaller. 
The results of Toivanen et al. differ substantially, as 
expected based on their poor reproduction of the spectra. 
They give 15%/55% smaller (S„) values for 129 Xe/ 131 Xe 
and both (S p ) values an order of magnitude smaller (for 
129 Xe also of opposite sign). We attribute these differ- 
ences to the sizeable truncations of the valence spaces in 
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FIG. 2: (Color online) Structure factors Soo, Sn and Soi for 
129 Xe as a function of u = p 2 b 2 /2, with harmonic-oscillator 
length b — 2.2853 fm. Results are shown at the one-body (lb) 
current level and including two-body (2b) currents with an 
estimate of the theoretical uncertainty given by the red band. 



FIG. 3: (Color online) Structure factors Soo, Sn an d Soi f° r 
131 Xe as a function of u = p 2 b 2 /2, with harmonic-oscillator 
length b = 2.2905 fm. The labels are as in Fig. H 



Refs. 1|| 21 1 and also because the nuclear interactions 
are not as well tested compared to this work. 

At finite p, one introduces isoscalar/isovector (0/1) 
structure factors Sbo(p), (p) and Sn(p) through 

Sa(p) = al Soo(p) + aoaiSoi(p) + a\ Su(p) . (17) 

Figures [2] and [3] show the predicted structure factors 
for 129 Xe and 131 Xe as a function of u = p 2 b 2 /2, with 
harmonic-oscillator length b = (h/mw) 1 ^ 2 and hw — 
(45A-V3 _ 25^ 2 / 3 ) MeV. Fits to S 00 , S n and S 01 are 
provided in the supplementary material. 

At the lb level at order Q 2 , the results shown corre- 
spond to the same current operator as in Refs. (l9l. l2lj|. 
but based on chiral EFT couplings and with the l/A 2 ^ 
term in gA and gp included. Since (S„) sets the scale 
for the p = response, one has Soo ~ Sii ~ — l/2«Sbi- 
Thus, our response is about 35% larger for 131 Xe com- 
pared to Ressell and Dean [21|, and about 30%/4 times 
larger for 129 Xe/ 131 Xe than Toivancn et al. 19]. For 
p > to,,-, u > 1.3, the structure factors nearly vanish for 
129 Xe, similar to Refs. [H, [2l|, as expected for a 1/2+ 
ground state, where only L = 1 multipoles contribute. 
For 131 Xe, the predicted structure factors remain nonva- 



nishing but small, similar to Ressell and Dean [21[. Both 
differ significantly from Toivanen et al. [19] suggesting 
very different higher L = 3 multipole contributions. 

When chiral 2b currents are included at order Q 3 , 
we provide theoretical error bands in Figs. [2] and [3] due 
to the nuclear uncertainties in WIMP currents in nu- 
clei. This takes into account the uncertainties in the 
low-energy couplings C3 , C4 (with the conservative range 
used in Ref. Q) as well as the range in nuclear densities 
p = 0.10. ..0.12 fm~ 3 . Chiral 2b currents provide impor- 
tant contributions to the structure factors, especially for 



p < 100 MeV, u < 0.7, where their effects are to signifi- 
cantly reduce, by about 25—55%, the isovector Sn factor, 
similar to the Gamow- Teller quenching Q and by about 
half of that the isoscalar-isovector Soi factor, compared 
to the results at the lb level at order Q 2 . As discussed, 
the 2b current contributions to the isoscalar Soo factor 
are small and are not considered here. The chiral 2b cur- 
rent contributions should be included in all calculations 
of spin-dependent elastic and inelastic WIMP scattering 
off nuclei. The effects of 2b currents are expected to be 
smaller for the vector structure factor and thus for the 
spin-independent elastic WIMP response. 

Summary. - This presents the first calculation of spin- 
dependent WIMP currents in nuclei based on chiral EFT, 
including the leading long-range 2b currents. They pre- 
dict a 25 — 55% reduction of the isovector part of the 
one-body axial-vector WIMP currents, where the range 
provides an estimate of the theoretical uncertainties of 
WIMP currents in nuclei. This should be included in 
limits on the WIMP couplings, where the spin-dependent 
analysis provides complementary constraints Q. As an 
application, we have calculated the structure factors for 
spin-dependent WIMP scattering off 129,131 Xe nuclei, us- 
ing the largest valence spaces accessible with nuclear in- 
teractions that have been tested in nuclear structure and 
double-beta decay studies in this region. Future work 
includes developing consistent interactions based on chi- 
ral EFT, where the present frontier is in the calcium re- 
gion 22], and investigating other nuclei and responses. 
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Supplementary Material 

TABLE I: Fits to the structure factors Soo, Su and Soi for spin-dependent WIMP elastic scattering off 129 Xe and 131 Xe nuclei, 
including lb and 2b currents as in Figs. [2] and [3] For the lb+2b current results, both the central value of the theoretical error 
band was used for the fits (first rows) and the limits of the band (second rows). The fitting function of the dimensionless 
variable u — p 2 b 2 /2 is an exponential e~ u multiplied by a ninth-order polynomial. The rows give the coefficients of the u n 
terms in the polynomial. 
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-0.122388 


-0.096693 


u 5 


—0 0133335 


-0.0179342 


-0.010598 


0.0317668 


0.0251835 


u 6 


00190579 

V/.V7V' 1 \J\> ' J • ./ 


0.00286368 


0.00168737 


-0.00492337 


-0.00392356 


u 7 


— 1 48373 x 10~ 4 


-2.65795 x 10~ 4 


-1.56768 x 10" 4 


4.39836 x 10" 4 


3.53343 x 10" 4 


u 8 


5.11732 x 10~ 6 


1.29656 x 10~ 5 


7.69202 x 10" 6 


-2.02852 x 10~ 5 


-1.65058 x 10 -5 


u 9 


-2.06597 x 10" 8 


-2.47418 x 10~ 7 


-1.48874 x 10" 7 


3.46755 x 10" 7 


2.88576 x 10" 7 


e~ u x 




Sn (lb+2b band) 


501 (lb+2b band) 


1 




0.0360513 


0.0226064 


-0.0888962 


-0.070431 


u 




-0.110705 


-0.0703558 


0.257562 


0.20644 


u 2 




0.150026 


0.0954932 


-0.336681 


-0.271749 


u 3 




-0.111714 


-0.0706386 


0.245328 


0.198772 


u 4 




0.0493115 


0.0308683 


-0.106746 


-0.0866783 


u 5 




-0.0130745 


-0.00810917 


0.0277603 


0.0226217 


u 6 




0.00208705 


0.00128522 


-0.00431334 


-0.00353705 


u 7 




-1.94174 x 10 -4 


-1.19086 x 10" 4 


3.86604 x 10" 4 


3.20472 x 10" 4 


u 8 




9.52295 x 10" 6 


5.84562 x 10~ 6 


-1.79013 x 10~ 5 


-1.51335 x 10~ 5 


u 9 




-1.83523 x 10 -7 


-1.13885 x 10 -7 


3.06893 x 10" 7 


2.70785 x 10" 7 








131 Xe 










u = p 2 b 2 /2, b = 2.2905 fm 




e~ u X 


Sm 


Su (lb) 


Su (lb+2b) 


Soi (lb) 


Soi (lb+2b) 


1 


0417889 


0.0368132 


0.022446 


-0.0784478 


-0.0608808 


u 


—0.111171 


-0.118361 


-0.0733931 


0.230484 


0.181473 


u 2 


171966 


0.176773 


0.110509 


-0.343106 


-0.272533 


U 3 


—0 133219 


-0.137987 


-0.0868752 


0.263525 


0.211776 


u 4 


0633805 


0.063821 


0.0405399 


-0.120946 


-0.0985956 


u 5 


— 0178388 


-0.0176743 


-0.0113544 


0.0331754 


0.027438 


u 6 


00282476 


0.00287653 


0.00187572 


-0.00528724 


-0.0044424 


u 7 


—2.31681 x 10~ 4 


-2.63605 x 10~ 4 


-1.75285 x 10" 4 


4.6475 x 10" 4 


3.97619 x 10" 4 


u 8 


7.78223 x 10~ 6 


1.23239 x 10" 5 


8.40043 x 10~ 6 


-2.00407 x 10 -5 


-1.74758 x 10~ 5 


u 9 


-4.49287 x 10" 10 


-2.17839 x 10~ 7 


-1.53632 x 10~ 7 


2.90375 x 10" 7 


2.55979 x 10" 7 


e -"x 




Su (lb+2b band) 


501 (lb+2b band) 


1 




0.0275713 


0.0173191 


-0.067911 


-0.0538485 


u 




-0.0891521 


-0.0576174 


0.20071 


0.162246 


u 2 




0.134094 


0.0868883 


-0.301325 


-0.24379 


u 3 




-0.105277 


-0.0684403 


0.233566 


0.190041 


u 4 




0.0490864 


0.0319771 


-0.108494 


-0.0887141 


u 5 




-0.0137279 


-0.00897632 


0.0301064 


0.0247694 


u 6 




0.00226103 


0.00148964 


-0.00485087 


-0.0040324 


u 7 




2.10252 x 10" 4 


-1.40246 x 10~ 4 


4.30628 x 10" 4 


3.64368 x 10~ 4 


u 8 




1.0004 x 10~ 5 


6.79325 x 10" 6 


-1.86527 x 10~ 5 


-1.62836 x 10~ 5 


u 9 




-1.80793 x 10" 7 


-1.26396 x 10" 7 


2.63842 x 10" 7 


2.48126 x 10" 7 



